Phase diagram of generalized fully frustrated XY model in two dimensions 



OO 

o 
o 

C 

a 
o 

(N 



Petter Minnhagen, 1 Beom Jun Kim, 2 Sebastian Bernhardsson, 1 and Gerardo Cristofano 3 

1 Dept. of Physics, Umea University, 901 87 Umed, Sweden 
2 Dept. of Physics, BK21 Physics Research Division and Institute of Basic Science, Sungkyunkwan Univ., Suwon 440-746, Korea 
3 Dipartimento di Scienze Fisiche, Universitd di Napoli "Federico II" and INFN, 
Sezione di Napoli, Via Cintia, Compl. universitario M. Sant'Angelo, 80126 Napoli, Italy 

It is shown that the phase diagram of the two-dimensional generalized fully-frustrated XY model 
on a square lattice contains a crossing of the chirality transition and the Kosterlitz-Thouless (KT) 
transition, as well as a stable phase characterized by a finite helicity modulus T and an unbroken 
chirality symmetry. The crossing point itself is consistent with a critical point without any jump 
in T, with the size (L) scaling T ~ £-°- 63 anc [ ^ e cr itical index v ~ 0.77. The KT transition line 
remains continuous beyond the crossing but eventually turns into a first-order line. The results are 
established using Monte-Carlo simulations of the staggered magnetization, helicity modulus, and 
the fourth-order helicity modulus. 
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I. INTRODUCTION 



The phase transitions of the two-dimensional (2D) 
fully-frustrated XY (FFXF) model on a square lattice 
has been a subject of controversy^' 3 -'^ The emerging 
consensus is that the model, as the temperature is low- 
ered, first undergoes an Ising-like transition associated 
with the chirality. At a slightly lower temperature it un- 
dergoes a universal jump Kosterlitz-Thouless (KT) tran- 
sition associated with the phase angles*^ Since the two 
transitions are extremely close to each other in tempera- 
ture the question whether there is only one merged tran- 
sition or actually two separate transitions has been the 
cause of the controversy. An argument for two separate 
transitions with the KT transition always at a lower tem- 
perature than the chirality transition was given by Kor- 
shunov in Ref. in terms of a kink-antikink instability of 
the domain walls separating domains with different chi- 
rality. The 2D generalized fully frustrated XY(GFFXY) 
model has the same degrees of freedom and the same 
symmetries as the FFXY model. The argument by Ko- 
rshuno\£ is quite general and appears to hinge only on 
the combined C/(l) and Zi symmetry of the model and 
the existence of Ising-like domain walls associated with 
the broken Z2 symmetry, 6 This strongly suggests that 
also the generalized model with the very same degrees of 
freedom and symmetry should always have a KT transi- 
tion at a lower temperature than the chirality transition. 
As shown here, this is not the case: The two transitions 
can merge in a single critical point. The reason for this 
unexpected result is that the symmetry of the model al- 
lows for a new phase 2D "quasi" phase-order. 

The present paper is organized as follows: In Sec. [Til 
we introduce the generalized fully-frustrated XY model. 
The results of our numerical simulations are presented in 
Sec. IIIII for the staggered magnetization, in Sec. IIVI for 
the helicity modulus, and in Sec. [V] for the fourth-order 
modulus, respectively. Finally, Sec. IVII is devoted to the 
summary of the paper. 



II. THE GENERALIZED FFXY MODEL 

The XY model on a square lattice in the presence of 
an external magnetic field transversal to the lattice plane 
is described by the action: 
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where di is the phase variable at the ith site, the sum is 
over nearest neighbors, J(> 0) is the coupling constant, 
T is the temperature, ks is the Boltzmann constant, 
and = (2e/hc) f. A • all is the line integral along the 
bond between adjacent sites i and j. We consider the 
case where the bond variables are fixed, uniformly 
quenched, out of equilibrium with the site variables and 
satisfy the condition Aij = 2irf: here the sum is over 
each set of bonds of an elementary plaquette and / is the 
strength of frustration. We assume that the local mag- 
netic field in Eq. (UJ is equal to the uniform applied field; 
such an approximation is more valid the smaller is the 
sample size L compared with the transverse penetration 
depth \±Z In the case of full frustration, i.e., / = 1/2, of 
interest to us here, such a model has a continuous U(l) 
symmetry associated with the rotation of spins and an 
extra discrete Zi symmetry, as it has been shown by an- 
alyzing the degeneracy of the ground stated Choosing 
the Landau gauge, such that vector potential vanishes on 
all horizontal bonds and on alternating vertical bounds, 
we get a lattice where each plaquette displays one antifer- 
romagnetic and three ferromagnetic bonds. Such a choice 
corresponds to switching the sign of the interaction. 

The generalized FFXY model is obtained by changing 
the form of the interaction from — JcoscA tc 10,11 
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This does not alter any symmetry present in the original 
FFXY model which corresponds to p = 1 since 2[1 — 



FIG. 1: The phase diagram of the 2D GFFXY model. The 
phases are characterized by the helicity modulus T and the 
staggered magnetization m. The four phases correspond to all 
possible combinations of finite and vanishing T and m. The 
points correspond to data obtained from the simulations. The 
horizontal line p = 1 corresponds to the usual FFXY model. 
The phase lines cross and merge in one point. 



cos 2 (0/2))] = 1 — coscj). The essential point is that U(<p) 
is periodic in 2ir and that the first term in an expansion 
for small <fi is second order, i.e., U(cf>) ~ J(j> 2 /2 for tfi <C 1. 

Figure [1] is the summary of the results present in this 
paper and shows the phase diagram in the (p, T)-plane 
as obtained from Monte-Carlo (MC) simulations. The 
helicity modulus T, which relates to the continuous an- 
gular symmetry, and the staggered magnetizations to, 
which relates to the discrete chirality symmetry, are used 
to detect phase boundaries. The phase diagram con- 
tains all four possible combinations of these two, i.e., 
(T, to) = (0, 0), (0, ^ 0), 0, 0), 0, ^ 0). The dashed 
horizontal line at p = 1 corresponds to the usual FFXY 
model, for which the phase (T ^ 0, m = 0) is not real- 
ized. 
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FIG. 2: MC determination of the staggered magnetization, 
(a) The transition across the horizontal phase line p c « 1.3479 
is confirmed by a sharp drop to zero for larger sizes L. (b) Size 
scaling of Binder cumulant B m for p = 1. The unique crossing 
of curves for different sizes yields the transition temperature 
T c h ~ 0.454. This method is used to determine the m-phase 
line for p < p c . 



III. STAGGERED MAGNETIZATION 

We first present numerical MC results of the staggered 
magnetization to, defined as 12 
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2=1 
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where (• • ■ ) is the ensemble average and the vorticity for 
the Ith elementary plaquette at (a;;, yi) is computed from 



si = (1/tt) 5Z{ij)ei 't'ij = w ith the sum taken in the 
anti-clockwise around the given plaquette. 

The ground states with the spontaneously broken chi- 
rality symmetry correspond to the two possible checker 
board patterns with alternating positive and negative 
vorticity. The energy per link in these ground states 
is given by U(ir/4) which corresponds to all links con- 
tributing the same energy. Since the two ground states 



with different checker board patterns are separated by 
an infinite energy barrier in the thermodynamic limit, 
the phase with the broken chirality symmetry persists at 
low enough temperatures as long as the pattern, where 
all links contribute the same energy, indeed corresponds 
to the ground state. However, this ceases to be true 
when p becomes larger than p c . In this new region the 
ground state instead corresponds to a pattern consisting 
of plaquettes with phase difference on three sides and 
7r on the remaining. The energy per link is hence instead 
U(ir)/A. The critical value p c is easily computed to be 
p c w 1.3479 from the condition that U(tv/4) = E7(tt)/4. 
Consequently, the new ground state at p > p c has no bro- 
ken chirality symmetry and hence corresponds to to = 0. 

Figure [Ha) illustrates the vanishing of the staggered 
magnetization to at p ps p c for T = 0.1 (the tempera- 
ture is in units of J/fcs throughout the present work). 
This horizontal part of the m-phase boundary eventually 
bends down toward smaller p as T is increased. This 
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FIG. 3: MC determination of the KT transition line, (a) The 
helicity modulus T as a function T for various sizes L. The 
broken lines correspond to the jump condition; b = 1 is the ex- 
pected universal jump. Tkt is estimated by extrapolating the 
crossing point with the 6 = 1 line and the data to L = oo, as 
shown in (b), which shows a second order polynomial extrap- 
olation of the crossing points and Tkt = 0.447 is obtained 
in accord with the finding Tkt = 0.446 in Ref. 0. (c) The 
difference AT = \TktQ>, L) — TktQj/2, L)\ in crossing points 
Tkt [see (a)] for (6,6/2) = (1,1/2) and (2,1). It is clearly 
seen that the jump has to be less than twice the universal 
jump consistent with the universal jump, (d) The same con- 
struction for p — 1.5. We believe that a nonunivcrsal jump 
for p > p c cannot be entirely ruled out. (a)-(c) correspond to 
p = 1 while (d) is for p = 1.5. 



part of the m-phase boundary we have traced out by the 
standard size scaling of Binder's cumulant B m for the 
order parameter m^ as displayed in Fig.[2Jb) for p = 1. 
T c h ~ 0.454 is obtained, in a good agreement with the 
earlier value 0.452 in Ref. 0. The complete m-phase line 
with marked data points are shown in Fig. [1] 



size L. These values are well approximated by a second 
order polynomial as shown in Fig. G2b). The extrapola- 
tion to L = oo gives Tkt — 0.447, which is very close 
to the value 0.446 obtained in Ref. 0. The close agree- 
ment shows that the method gives a good estimate of the 
KT transition temperature. The data points in Fig.[T]for 
p < 1.32 are obtained by this method. One notes that 
the m-phase line and the KT-line are extremely close for 
these p- values and only the smaller p- values, like p = 0.5, 
display a clear separation within our accuracy. 

The determination of the KT-line rests on the assump- 
tion that the KT-jump has the universal value 2/ir. A 
jump means that the crossing point between the T(T) 
and b(2/ir)T should give the same Tkt(L — oo) for all 
b < 1. Figure 02c) shows the difference AT(b,L) — 
TKr(b,L) — TKT(b/2, L) for 6=1. Our result is consis- 
tent with a universal jump KT transition, since AT(6 = 
l,L) is consistent with a vanishing for L — oo. On the 
other hand the jumps size is inconsistent with a double 
jump since AT(b — 2, L) approaches a finite value. For 
larger values of p > p c , like p = 1.5, the jump is, on 
the other hand, consistent with a jump larger than the 
universal jump as illustrated in Fig. [3fd): For this value 
AT(6 = 2, L) is consistent with a vanishing, suggest- 
ing that the jump at the KT transition could be larger 
than the universal KT value. The transition at p = 1.5 
shows no sign of any first order character from which we 
conclude that it is continuous. In this case the jump is 
expected to have the universal value. Our data neither 
support nor rule out this expectation. Conversely, a con- 
tinuous KT transition with a nonuniversal jump can nei- 
ther be ruled out. However, when p is increased further 
the transition does eventually become first order as can 
be detected from the double well structure in the energy 
histogram. For the first order transition the jump should 
be nonuniversal and larger than the universal jumpiiS 
In the limit of p = oo the model reduces to the infinite 
state Potts models which is known to have a first order 
transition. 



FOURTH-ORDER MODULUS 



IV. HELICITY MODULUS 

The quasi 2D phase ordering is measured by the he- 
licity modulus T defined as the stiffness in response to 
the twist 5 of the phase variables across the system: 
T = (d 2 F/dS 2 )s=o where F is the free energy. The condi- 
tion for a KT transition is characterized by the universal 
jump in the helicity modulus, T(Tkt)/Tkt = 2/7T i 14 i 17 . 
Thus a KT transition can be located by the crossing point 
between the line y — (2/tt)T and the helicity modulus 
curve y = T(T) as illustrated for p = 1 in Fig. [3fa). 
In practice, a precise determination requires the difficult 
task of extrapolating to L = ooi Here we use the follow- 
ing method: The values of the Tkt at the crossing point 
with the line y = (2/tt)T are determined as a function of 



According to the argument given by Korshuno\«* the 
KT transition always occurs at a lower temperature than 
the chirality transition. This is consistent with the result 
we find for p- values below the horizontal line in the phase 
diagram. In this part of the phase diagram the argument 
by Korshunov a is valid and the chirality transition and 
the KT transitions are separated with the KT transi- 
tion always at a lower temperature. On the other hand, 
when the horizontal phase line meets and crosses the KT 
line, Korshunov's argument^ is no longer valid and one 
expects a merged character of the transition. The ar- 
gument by Korshunov fails because it presumes the ex- 
istence of Ising-like domain walls and such walls do not 
exist above the horizontal line. In order to monitor the 
change of character of the transition we study the fourth- 
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FIG. 4: MC simulation of the fourth-order modulus T4.(a) 
shows the typical T4-characteristic of a KT-transition for 
p — 1: The well minimum moves from higher T towards the 
transition temperature (vertical broken line) with increasing 
L. The size of the minimum extrapolates to a finite value, (a) 
and (b) show that this typical KT-feature remains intact as 
p is increases to the vicinity of p c . (c) and (d) show that this 
KT-feature is dramatically changed in the immediate vicinity 
of p c , while (e) and (f) show that the KT-feature is recovered 
for p- values above p c 



order helicity modulus, 15 defined by the expansion of the 
free energy AF = F(S) - F(0) = TS 2 /2! + T 4 <5 4 /4!. The 
observation that T4 is finite and negative precisely at the 
KT transition 15 leads to the conclusion that T makes an 
abrupt jump at Tkt i n order to fulfill the requirement 
of AF > at any T. Thus T4 offers a way to verify a 
discontinuous jump at a KT transition^ 

Figure[4ja) for p = 1 is consistent with the typical KT- 
features for T4: The minimum position well approaches 
Tkt from above (the vertical lines in Fig. 2] mark the ex- 
pected positions of the transition) and the depth of the 
well remains finite, as L is increased, confirming the ex- 
istence of an abrupt jump of the helicity modulus?^ The 
shift of the minimum position towards lower tempera- 
tures constitutes the characteristics of a KT transition 
and is well established up to p — 1.31 [see Fig. B|b)]. 
As p is increased through the crossing region around 
p c w 1.3479 there are dramatic changes but for larger p 
the typical KT behavior of T4 reappears [see Fig. H£e)]. 
This is consistent with a crossing where a KT transition 
disappears and reappears as p is increased. The charac- 
teristics close to p c ss 1.3479 is instead consistent with 
T4 = as L — > 00 [see Fig. |4jc) and (d)] . 

If T4 = then the helicity modulus T does not need 
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FIG. 5: MC-simulation of critical size scaling for the helicity 
modulus T at the crossing point of the phase lines, (a) The 
size scaling T ~ L~ a is obeyed to very good approximation. 
The collapsing point of the data determines the critical tem- 
perature T c » 0.1675. (b) With the T c from (a) the fuller 
scaling T = L~ a F[L 1/v (T - T c )] with v = 0.77 is born out. 



to have a jump at the transition, which opens up the 
possibility of a continuous vanishing of T and the critical 
scaling T ~ L~ a . FigureG^a) shows that such a size scal- 
ing is indeed obtained close to p c ~ 1.3479. Furthermore, 
Fig. E^b) shows that the standard critical scaling form for 
a continuous phase transition T = L~ a F\L 1 l v '{T — T c )] 
is also valid to very good approximation which suggests 
that the correlation length £ diverges as £ ~ \T — T c \~ v ■ 
The values obtained for the critical indices are a ss 0.63 
and v w 0.77. 



VI. SUMMARY 

The main result of the present work is that in general 
a model with the same symmetry and degrees of free- 
dom as the 2D fully frustrated XY model can have four 
stable phases. Only three of these phases are present in 
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the usual FFXY model. The new phase allowed by sym- 
metry combines unbroken chirality with quasi 2D phase 
ordering. The existence of this phase also means that the 
KT and chirality phase lines cross. The crossing point is 
a critical point at which the helicity modulus obeys scal- 
ing and vanishes smoothly without a universal jump. 
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